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We show that, in contrast to known results in the massive case, a minimally gauged mass¬ 
less Rarita-Schwinger field yields a consistent classical theory, with a generalized fermionic 
gauge invariance realized as a canonical transformation. To simplify the algebra, we study a 
two-component left chiral reduction of the massless theory. We formulate the classical theory 
in both Lagrangian and Hamiltonian form for a general non-Abelian gauging, and analyze 
the constraints and the Rarita-Schwinger gauge invariance of the action. An explicit wave 
front calculation for Abelian gauge fields shows that wave-like modes do not propagate with 
superluminal velocities. An analysis of Rarita-Schwinger spinor scattering from gauge fields 
shows that adiabatic decoupling fails in the limit of zero gauge field amplitude, invalidating 
various “no-go” theorems based on “on-shell” methods that claim to show the impossibil¬ 
ity of gauging Rarita-Schwinger fields. Quantization of Rarita-Schwinger fields, using many 
formulas from this paper, is taken up in the following paper. 


I. INTRODUCTION 


A. Motivations and Background 


Cancelation of gauge anomalies is a basic requirement for constructing grand unified models, 
and the usual assumption is that anomalies must cancel among spin ^ fermion fields. However, 
a 1985 paper of Marcus [1] showed that in principle an SU (8) gauge theory can be constructed 
with spin | Rarita-Schwinger fermions playing a role in anomaly cancelation, and we have recently 
constructed [2] a family unification model incorporating this observation. Using gauged spin | 
fields in a grand unification model raises the question of whether such fields admit a consistent 
quantum, or even classical theory. It is well known, from papers of Johnson and Sudarshan Q] 
and Velo and Zwanziger jd], and much subsequent literature (see e.g. Hortacsu [5] and Deser and 
Waldron jf|), that theories of massive gauged Rarita-Schwinger fields have serious problems. Does 
setting the fermion mass to zero eliminate these difficulties? 

The lesson we have learned from the success of the Standard Model is that fundamental fermion 
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masses lead to problems and are to be avoided; all mass is generated by spontaneous symmetry 
breaking, either through coupling to the Higgs boson or through the formation of chiral symmetry 
breaking fermion condensates. So from a modern point of view, the Rarita-Schwinger theory 
with an explicit mass term is suspect. Several hints that the behavior of the massless theory 
may be satisfactory are already apparent from a study of the zero mass limit of formulas in the 
Velo-Zwanziger paper. First, in their demonstration of superluminal signaling, the problematic 
sign change that they find for large B fields (Eq. (2.15) of 4]) is not present when the mass 
is set to zero. Second, when the mass is zero, the secondary constraint that they derive (Eq. 
(2.10) of 4|]) appears as a factor in the change in the action under a Rarita-Schwinger gauging 
5tpfj, = with D^ the usual gauge covariant derivative. (This statement is not in [4], but is 
an easy calculation from their Eqs. (2.1)—(2.3), with the D^ of this paper their Hence, 

the constrained action in the massless gauged Rarita-Schwinger theory has a fermionic gauge 
invariance that is the natural generalization of the fermionic gauge invariance of the massless free 


Rarita-Schwinger theory. Third, their formula for the anticommutator (Eq. (4.12) of So in the 
zero mass case develops an apparent singularity in the limit of vanishing gauge field B, and so 
their quantization does not limit to the standard free theory quantization. However, since the 
massive theory does not have a fermionic gauge invariance, Ref. 4] does not include a gauge-fixing 
term analogous to that used in the massless case, but gauge fixing is needed to get a consistent 
quantum theory for a free massless Rarita-Schwinger field. So these observations, following from 
the equations in 4[, suggest that a study of the massless Rarita-Schwinger held coupled to spin-1 
gauge fields is in order. 

In a different and more recent setting, massless Rarita-Schwinger fields appear consistently 
coupled to gravity as the gravitinos of supergravity, as discussed by Das and Freedman [7]. Grisaru, 
Pendleton and van Nieuwenhuizen [8] have shown that soft spin | fermions must be coupled to 
gravity as in supergravity, in an analysis based on the free particle external line pole structure of 
spin | fields that do not have spin 1 gauge couplings. Their result has been extended to gauged 


spin | fields in various recent “no-go” theorems based on “on-shell” methods [9], Q, that again 
assume a free particle external line pole structure. None of these papers have analyzed the gauged 
Rarita-Schwinger equation to determine the asymptotic held structure. Thus, these papers do not 
prove that there cannot be a consistent theory of massless, gauged Rarita-Schwinger helds, so again 
a detailed study of this possibility is warranted. 
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B. Outline of the paper, and summary 


With these motivations and background in mind, we embark in this paper on a detailed study of 
the classical theory of a minimally gauged massless Rarita-Schwinger field. In Sec. 2, we give the 
Lorentz covariant Lagrangian for a gauged four-component Rarita-Schwinger spinor field, derive 
the source current for the gauge field, and check that it is gauge-covariantly conserved. We also 
give the Lorentz covariant form of the constraints, of the fermionic gauge transformation, and 
of the symmetric stress-energy tensor, and briefly discuss the generalization to non-flat metrics. 
Since in the massless case left chiral and right chiral components of the field decouple, in Sec. 3 
we rewrite the Lagrangian for left chiral components in terms of two-component spinors and Pauli 
matrices, which simplifies the subsequent analysis. We then give the Euler-Lagrange equations 
in two-component form, and use them to analyze the structure of constraints and the fermionic 
gauge transformation of the action. In Sec. 4 we introduce canonical momenta for the Rarita- 
Schwinger field components, which are used to define classical Poisson brackets, and discuss the 
role of the constraints as generators of gauge transformations under the bracket operation. We 
show that the constraints group into two sets of four, within each of which there are vanishing 
Poisson brackets. In Sec. 5, we argue that fermionic gauge transformations give a generalized 
form of gauge invariance, corresponding to the presence of redundant gauge degrees of freedom, by 
studying the properties of both infinitesimal and general finite gauge transformations. We show 
that infinitesimal gauge transformations are an invariance of the constrained action functional that 
governs the influence of Rarita-Schwinger fields on gauge and gravitational fields. We show that 
finite gauge transformations take the form of generalized auxiliary fields, which lead to an extended 
action that has an exact invariance under fermionic gauge transformations. In Sec. 6 we specialize 
to the case of an Abelian gauge held (as in Q) and analyze the wavefront structure, showing that 
physical wave modes propagate with luminal velocities; an extension of this discussion, showing 
that gauge modes are subluminal is given in Appendix B. In Sec. 7, making a transition to first 
quantization, we analyze Rarita-Schwinger fermion scattering from an Abelian gauge potential. 
We show that the asymptotic state structure assumed in “on-shell no-go” theorems is not realized, 
but that a consistent scattering amplitude can be formulated using an analog of the distorted wave 
Born approximation. In Sec. 8 we give a brief summary and discussion, and in Appendix A we 
summarize our notational conventions and some useful identities. We suggest that the reader skim 
through Appendix A before going on to Sec. 2, since things stated in Appendix A are not repeated 
in the body of the paper. In the paper that follows this one we build on our analysis to discuss 
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quantized Rarita-Schwinger fields. 

II. LAGRANGIAN AND COVARIANT CURRENT CONSERVATION IN 

FOUR-COMPONENT FORM 

A. Flat spacetime 

The action for the massless Rarita-Schwinger theory is 

SM =\I , 

R mu =ie^ p (^Xp(D^ p ) u , 

=d v ^ + , 

A U _ A A.U 
■^■vv u Av 5 

(i) 

with ^ au = 'ip liau (x, t) a four-vector four-component spinor, with four-vector index g = 0, 
spinor index a = 1,4, and SU(n ) internal symmetry index u = with SU(n ) gauge 

generators tA, A = l,...,n 2 — 1. Taking u to range from 1 to n means that, for definiteness, 
we are assuming that the spinors transform according to the fundamental representation of the 
SU(n ) internal symmetry group, but other representations and other compact Lie groups can be 
accommodated by assigning the internal indices u and A the appropriate range. Note that tA, A™ v , 
and D u all commute with the gamma matrices and the Pauli spin matrices from which the gamma 
matrices are constructed, and for an Abelian internal symmetry group, the indices u and A are 
not needed. Using 

VW = Vi/juifrVa . ( 2 ) 

together with the adjoint convention (x\x 2 )^ = X 2 X 1 f° r Grassmann variables yi, X 2 , it is easy to 
verify that S is self-adjoint. 

From here on we will usually not indicate the spinor indices a , /3 and internal symmetry indices 
u, v explicitly, but they are implicit in all formulas. Varying S with respect to the Rarita-Schwinger 
fields, we get the equations of motion 

^ vp dA plr] =ge^Xt A x v , 

( 3 ) 
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Re-expressed in terms of the covariant derivative, these are 

7 „=0 , 

t mvp lr] D v ^ p = 0 . 

(4) 


The g = 0 component of these equations gives the primary constraints 

7e =0 , 
e enr le D n ^ r = 0 , 

(5) 

with e,n,r summed from 1 to 3. Contracting the equation of motion for with g^^D p and the 
equation of motion for ip p with g~ 1 D p , we get the secondary constraints 

e^ pFpu7v = 0 , 

e mup l v F^ p =0 , 

( 6 ) 


where we have introduced the gauge held strength 

F pu =g~ 1 [D p ,D 1/ \ = g^lSp,^5 u \ 
= d p A u — d v A p + g[A p , A v \ , 


(7) 


which with the adjoint representation index A indicated explicitly reads 


B aC 


F* = d p A* - d u A* + gfABcA^A 


( 8 ) 


Under a Rarita-Schwinger gauge transformation (with e a four-component spinor), which is a 
natural gauge held generalization of the fermionic gauge invariance for a free, massless Rarita- 


Schwinger held discussed in 




Vv , fcVv = F> p e , 

Vv , 


(9) 
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the action of Eq. © changes according to 


= -\i g J dA x + 0{e...e) . (10) 


The factors bracketed in large parentheses are identical to the secondary constraints of Eq. ([6]). 
This equation holds with finite (not necessarily infinitesimal) e and its adjoint eq the precise form 
of the quadratic term is given in Eq. (|72j) below. We will argue in Sec. 5 that Eq. (HOD implies 
that, even when coupled to gauge fields, the Rarita-Schwinger theory has a generalized form of 
fermionic gauge invariance. 

Adding the gauge field action 


S(A A ) = ~j d 4 xF A F A ^ , (11) 

and varying the sum Styfj,) + S(A A ) with respect to the gauge potential, we get the gauge field 
equation of motion 


D v F a » v = d v F A ^ + gfABcA^F 0 ^ = gJ A » , 

J Afl rjtA'I’p ■ 

( 12 ) 


A straightforward calculation using Eqs. ([3|) shows that the gauge field source current J Afl obeys 
the covariant conservation equation 


+ gfABcAp c » = 0 


(13) 


as required for consistency of Eq. (1121) . So from the Rarita-Schwinger and gauge field actions, we 
have obtained a formally consistent set of equations of motion. 

In addition to the gauge field source current, there is an additional current J^ that obeys an 
ordinary conservation equation, 


jm 

duJ 11 =0 . 


(14) 


In the massive spinor case, Velo and Zwanziger 


4|] argue that the analogous current, within the 


constraint subspace of Eq. (j5|) , should have a positive time component. In the massless case 
we see no reason for this requirement, since Eq. (1141) is the fermion number current, and its time 
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component, giving the fermion number density, can have either sign. However, we shall use parts of 
the positivity argument of 4] later on in discussing positivity of the Dirac bracket anticommutator. 

The symmetric stress-energy tensor for the free massless Rarita-Schwinger has been computed 
by Das [12] (see also Allcock and Hall {3]). Changing ordinary derivatives to gauge covariant 
derivatives, Das’s formula becomes 


rpCTT _ ^ bXllVO 

J RS — - 


^A75(7X +l a K) D ^P 




+^d a yl’\ibi^[i a + b a -,i T WM P 


(15) 


(This formula can be made manifestly self-adjoint by replacing D u by \{D U — D u ), but this is not 
needed to verify stress-energy tensor conservation.) Adding the gauge field stress-energy tensor, 

Tgbuge = ~\v aT F^F AX » + F a "F ax t , (16) 

a lengthy calculation, using Eq. (1131) together with identities and alternative forms of the equations 
of motion given in Appendix A, shows that the total tensor is conserved, 


d ff {T&+Tr 


gauge. 


= 0 


(17) 


B. Generalization to general gp V 

The generalization of the Rarita-Schwinger action to curved spacetime has been reviewed by 
Deser and Waldron [6]. In Eq. ([!]), d^x is replaced by the invariant volume element d 4 x(— g) 1 ^ 2 , 
and the covariant derivative D u becomes the curved spacetime covariant derivative 

D v i\) p = d v i\) p - r^/3 + + gA v ^ p , (18) 

with T„p and u; umn the affine and spin connections. The Rarita-Schwinger equation of Eq. (HD and 
the primary constraint of Eq. (JSD have the same form as before, in terms of the extended covariant 
derivative D u . The secondary constraint of Eq. Q now reads 

e^V p [^,K]7r,=0 , 

e^ lr] [Dp,D^ p = 0 , 


(19) 
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with u v defined by the adjoint of D u . The commutator of covariant derivatives is now given by 

fl 


[Dp, D u ]^ p = -R^pl/Ja + jRpumnl^^p + Q^pv^p 


( 20 ) 


with RJ lup and Rpvmn components of the Riemann curvature tensor, and as in flat spacetime involves 
only ijjp and not its time or space derivatives. In terms of the extended covariant derivative, the 
fermionic gauge transformation is still given by Eq. Q, and under this gauge transformation the 
change in the action is now given by 


= ~i J d A x e7 5 (e w ^[I} M , A,]Vv) + (f nvp 'P P ^n^vYlv)^ +0(e...e) , (21) 


with the factors bracketed in large parentheses now identical to the secondary constraints of Eq. 
(1191) (and again with e and finite). The arguments to be given in Sec. 5 then imply that in the 
presence of both gravitation and gauge fields, the Rarita-Schwinger theory has a generalized form 
of fermionic gauge invariance. Having established this curved spacetime generalization, we will 
continue in the remainder of this and the following paper to work in flat spacetime, but we expect 
everything done in what follows to have a curved spacetime generalization when the covariant 
derivative is suitably extended. 


III. LAGRANGIAN ANALYSIS FOR LEFT CHIRAL SPINORS IN 
TWO-COMPONENT FORM 

Although we could continue with the four-component formalism to study constraints, the Hamil¬ 
tonian formalism, and quantization, it will be more convenient to first reduce the four component 
equation to decoupled equations for left and right chiral components of 0“ (with a the spinor 
index and with the internal symmetry index implicit). Since these are related by symmetry, we 
can then focus our analysis on the two-component equations for the left chiral component, which 
is the component conventionally used in formulating grand unified models (see, e.g. Q]). 

We convert the action of Eq. © to two-component form for the left chiral components of 0“, 
using the Dirac matrices given in Eqs. (1A2I) and (1A4D . Defining the two-component four vector 
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spinor and its adjoint by 

PrM^ = ( ^ m = 0,1,2,3 , a = 1,2 , 

^la F L=(^la 0 ) , 


( 22 ) 


the action decomposes into uncoupled left and right chiral parts. The left chiral part, with spinor 
indices a suppressed, is given by 

S(^) = i J d 4 x[-^ 0 a ■ D x T + • a x Z)To + • D x T - • a x £> 0 tf] . (23) 

Varying with respect to we get the Euler-Lagrange equation 

0 = V = a x DTo + .D x ^ - a x D 0 ^ , (24) 

while varying with respect to Tq we get the primary constraint (given in four-component form in 
Eq. ©) 

0 = Vo = x = ^- Dx^ . (25) 

(The abbreviation Vq = x conforms to the notation of [d|.) A second primary constraint follows 
from the fact that the action has no dependence on dTg/df, which implies that the momentum 
conjugate to Tq vanishes identically, 


P ^=° 


(26) 


Contracting V with a and with g l D, and using the covariant derivative relations of Eq. (IA14I) . 
we get respectively 


<7 • V =2 id + y , 

g^D -V =iu + g~ l D 0 x , 


(27) 


with 


6 =<j ■ D'&q — Dqg • T , 
oj =a ■ o — (B + a x E) • T 


( 28 ) 
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Since the Euler-Lagrange equations imply that V and x vanish for all times, we learn that 6 and uo 
vanish also for all times. Since 9 involves a time derivative, its vanishing is just one component of 
the equation of motion for . But io involves no time derivatives, so it is a secondary constraint 
that relates d'o to (given in four-component form in Eq. ([6])). For each of the above equations, 
there is a corresponding relation for the adjoint quantity. 

The equation of motion V = 0 can be written in a simpler form by using the identities of Eqs. 
(jAlOl) and (IA11I) as follows. Using Eq. ()A10f) to simplify 0 = <7 x V — iV, we get an equation for 

A)*, 

D 0 $ = D\ h 0 + ~ [— a x (D x T) + ID x \k] . (29) 

A further simplihcation can be achieved by incorporating the primary constraint x = 0, through 
applying Eq. (1A11D to A = D x T, 

0 = ax = vv-(.Dx$) = Dx$-iax(Dx$) . (30) 

Using this to replace the first term in square brackets in Eq. (|29l) we get the alternative form of 
the equation of motion, valid when the constraint x = 0 is satisfied, 

£) 0 T = D^ 0 + iD xl . (31) 

Writing the gauge field interaction terms in Eq. (1231) in the form 

Sint(’JV) = f / d*x(AZj m + A b • J B ) , (32) 

we find the left chiral contribution to the current of Eq. (1121) in the form 

jAO = _ . g x ^ s 

J A =^ 0 tAd x T + x at a 'ho - x tA$ ■ 

(33) 


Replacing t,A by —i, we find the corresponding singlet current in the form 

J° =*Tt • a x T , 

J = - i(Tj<r xl + l*x <tT 0 - x f ) . 

(34) 


For the energy integral computed from the left chiral part of the the stress-energy tensor of Eq. 
(fl5l) . we find 


d 3 xT™ = - 


1 


H = - 


• D x T 


(35) 
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To conclude this section, we verify that the action of Eq. (1231) has a fermionic gauge invariance 
on the constraint surface u = 0 , = 0, as already seen in covariant form following Eq. ([9]). 

Letting e be a general space and time dependent two-component spinor, we introduce the fermionic 
gauge changes 

T —+ <5 G T , <5 G T = De , 

*0 —+ ‘Sg'I'o , 5 g T 0 = D 0 e , 

(36) 

and their adjoints, which are the left chiral form of the gauge change of Eq. Q. Substituting this 
into Eq. (f23l) . integrating by parts where needed, and using Eqs. (IA14I) to simplify commutators 
of covariant derivatives, we find that Eq. (1101) takes the two-component spinor form 

ScSi^n) = J d 4 x(uj*e - eV) + 0(J...e) , (37) 

with the quadratic term given in Eq. (1701) below. Hence the action on the constraint surface 

cj = = 0 has a fermionic gauge invariance. Another gauge invariant, on the constraint surface 

X = X 1 = 0) i s the fermion number, given by the space integral of the time component of the singlet 

current of Eq. (|34D . f d 3 xJ °, which has the gauge variation 

5g j d 3 xJ° = J d 3 x[—i(e <i x + y^e) + ge^d ■ Be] . (38) 

Again, these equations hold for e and it adjoint finite. 

However, neither the equation of motion, the constraints x and w, the non-Abelian “charge” 
f d 3 xJ B0 , nor the integrated Hamiltonian H are gauge invariant in the interacting case. Using <5 G 
to denote gauge variations, we have 

5qV = — ig(B + (j x E)e , 

5g0 = — iga ■ Ee , 

$GX = -igd ■ Be , 

<5 g w =a ■ BD^e — (B + (T x E) ■ De , 

Sg J d 3 xJ B() =g J d 3 x(^[A, ts\ ■ a x T + x a ■ [Lb, A]e^ , 

6 g H =\ig f d 3 x{^ ] • Be - e ] B ■ T) . 


(39) 
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The only global fermionic gauge invariants are the action integral, and the fermion number integral, 
in both flat and curved spacetimes. 

These results have an interpretation in terms of the distinction between a gauge transforma¬ 
tion, customarily defined as an invariance of the physical state of the system, and a canonical 
transformation. The usual gauge transformations in gauge field theories and general relativity are 
invariances of the action without imposition of a constraint, and consequently are invariances of 
the field equations and the Hamiltonian. Such gauge transformations are a special case of canon¬ 
ical transformations, but the converse is not true: canonical transformations in general alter the 
action, the field equations, and the Hamiltonian. We will see in Sec. 4 that the fermionic gauge 
transformation of Eq. (1361) are always canonical transformations, which reduce to gauge transfor¬ 
mations of the customary type only when the external gauge fields vanish. However, by virtue of 
the Jacobi identity for the Poisson bracket, canonical transformations preserve inner properties of 
the theory. As an example, that will be needed in our further discussion of generalized fermionic 
gauge invariance in Sec. 5, we verify that the secondary constraint following from the gauge-varied 
equation of motion V and primary constraint Vq = X agrees with the gauge variation of the original 
secondary constraint u. From Eq. CUD we have 

D ■ V - D 0 x = igu . (40) 

Preservation of inner properties under the fermionic gauge transformation means that we should 
find that 

D ■ 5 g V - D 0 5 g X = ig^GU • (41) 

Substituting Eqs. (1391) into the left hand side of Eq. (1411) gives 

ig[Do<j ■ Be — D ■ (B + a x E)e)\ = ig[a ■ BDqc — (B + a x E) ■ De + Ce] , (42) 

with the commutator remainder C given by 

C = a- [D 0 B - BD 0 + DxE + ExD]-(D-B-B-D) = 0 , (43) 

which vanishes by virtue of the gauge field Bianchi identity. 

In Sec. 5 we will discuss in more detail why the fermionic gauge transformation, because it leaves 
the constrained action invariant, corresponds to an unwanted redundancy in the time evolution. 
To break the gauge invariance we can introduce an additional constraint, in the form 


/W = 0 


( 44 ) 
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with / a scalar function of its argument. This constraint, together with the x constraint, leaves one 
independent two-component spinor of the original three in T, corresponding to the physical massless 
Rarita-Schwinger modes propagating in the gauge field background. We will limit ourselves to 
considering linear constraints of the general form 

f = L-$ , (45) 

and the choice L = D, a gauge covariant radiation gauge analog, plays a special role in our analysis. 
By not specializing L in our formulas, we can also examine the consequences of omitting a gauge 
fixing condition, corresponding to taking L = 0. 

We proceed to examine the gauge covariant radiation gauge condition in more detail. We note 
that since 

a ■ Da ■ = D ■ + ix, (46) 

the primary constraint x = 0 implies that 

a ■ Da • # = D • # . (47) 

Hence when a ■ D is invertible, which is expected in a perturbation expansion in the gauge coupling 
g, the covariant radiation gauge constraint D ■ 4/ = 0 implies that 

a • T = 0 . (48) 

Conversely, Eqs. (1461) and (1471) show that D • T = 0 and a • T = 0 together imply the primary 
constraint x = 0, and also a ■ Ik = 0 and X = 0 imply D ■ ik = 0. 

We next note that on a given initial time slice, covariant radiation gauge is attainable. Under 
the gauge transformation of Eq. (l36l) , we see that 

D ■ T -> D ■ T + {Dfe . (49) 

Hence when (D) 2 is invertible, which we expect to be true in a perturbative sense, then we can 
invert (D) 2 e = —D ■ T, to find a gauge function e that brings a general T to covariant radiation 
gauge. Since 

{a-Df = (Df+ga-B , (50) 

the conditions for a ■ D to be invertible, and for ( D) 2 to be invertible, are related. For generic 
non-Abelian gauge fields both of these operators should be invertible, but there will be isolated 
gauge field configurations for which a ■ D has zeros. 
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However, although covariant radiation gauge can be imposed on any time slice, it is not preserved 
by the equation of motion for T. To see this, let us consider the simplified case in which the gauge 
potential is specialized to = 0 and 8qA = 0, so that only a static B field is present. From Eq. 
m we have 

d 0 (D-$) = {D) 2 * 0 + gB-$ = [(D) 2 + ga-B]* 0 = (&-D) 2 * 0 . (51) 


So do (D ■ T) = 0 implies To = 0, but this is one constraint too many. Hence at each infinitesimal 
time step, we must make a further infinitesimal fermionic gauge transformation to maintain the 
covariant radiation gauge condition, as further discussed in Sec. 5B below. Only in the absence of 
gauge fields can we simultaneously impose the constraints V • ip = 0, a ■ ip = 0, and ip o = 0, as used 
in the discussion of 11] for the free Rarita-Schwinger case, 


IV. CANONICAL MOMENTA, CLASSICAL BRACKETS, AND GAUGE 

GENERATORS 


We next introduce the canonical momentum conjugate to T, defined by 

B dLs 1 ft - 
P =-— = — T' x a , 

9(9 0 T) 2 

which can be solved for 4H using the final line of Eq. (IA11D . 


T f = iP - P 


X (T 


(52) 


(53) 


We will use Eq. (15311 when computing classical brackets involving \fd using the formula of Eq. 
(1A17I) . Eq. (1521) can be written as an explicit matrix relation for the six components of P and T^, 


(Pi) 


( 0 

0 

1 

0 

0 



( 

pi 


0 

0 

0 

-1 

i 

0 



P r 

■*2 

1 

-1 

0 

0 

0 

0 

1 



pi 

“ 2 

0 

1 

0 

0 

1 

0 



P t 


0 

i 

0 

-1 

0 

0 



VdJ 


v -* 

0 

-1 

0 

0 

0 ) 




( 54 ) 


showing that they are related by an anti-self-adjoint matrix with determinant —1/16. 
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The four constraints introduced in Sec. 3 are 


=P *l ’ 

02 =(<? • S) _1 w = \ho — (<? ■ -B) _1 (-B + ex x E) ■ , 

03 =X = v ■ D x \k , 

0 4 =L ■ T 

(55) 


In writing these we are assuming that a ■ B is invertible in the non-Abelian case. We are writing 
the gauge fixing condition as a general linear gauge fixing constraint L ■ T, so as to keep track of 
which terms in the final answers arise from gauge fixing, which is not evident if we specialize by 
replacing L by D at this stage. The constraints of Eq. (1551) . including the gauge fixing constraint 
04, are all first class in the Dirac classification, since they have vanishing mutual classical brackets. 
This is a consequence of the fact that starting with a constraint depending on T but not on 
and taking an arbitrary number of time derivatives, one still has a constraint depending only on 
T. 

To preserve the adjoint properties of the Rarita-Schwinger equation, for each of these four 
constraints we must impose a corresponding adjoint constraint. Using Eq. (|53l) to express 'Id in 
terms of P, we write these as 


Xi =C p *t) t = ~ P *o i 

X 2 =w t (a-R)“ 1 = Tj -P- [i(B + a x E) - & x (B + a x E)\(a ■ B)~ l , 

X 3 =X f = 2 P-t) , 

X 4 ■ L = P ■ {iL — a x L) 

(56) 


(The reason for the minus sign in the definition P ^t = — P^, Q will be given in Sec. 2 of the following 
paper where we discuss the Hamiltonian form of the equations.) The constraints 0 a are implicitly 
2 n component column vectors, and the adjoint constraints Xa are implicitly 2 n component row 
vectors, with 2 n arising from the product of a factor of 2 for the two implicit spinor indices, and a 
factor of n for the n implicit SU (n) internal symmetry indices. 

When L = D, we see that 0 4 becomes 04 = D ■ \k, and xa becomes X 4 = iP - ^D — P • & /Id = 
(*/ 2 )x 3 — P ■ a x *D. So a special feature of covariant radiation gauge, which will be exploited later, 
is that the constraints 03, 0 4 are contractions ofifxD and D with T, and the constraints X 3 , X 4 
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are contractions of linear combinations of the duals D and a x D with P. That is, in covariant 
radiation gauge the constraint spaces selected by X3>X4 and ^> 3,04 are duals of one another. 

We can now compute the classical brackets of the constraints. We see that the brackets of the 
4>s and x s vanish among themselves, 


\4>aAb\c =0 * 

[Xa,Xb\c = 0 , 
a, b = 1 ,..., 4 


( 57 ) 


On the other hand, the brackets of the cj)s with the x s give a nontrivial matrix of brackets M, 
which has a nonvanishing determinant, 


M a b{x, y) =[4> a (x),Xb{y)]c + 0 , 
det M . 


( 58 ) 


Thus, in terms of the Dirac classification, the original first class constraints (j) a have become second 
class, not from adding new constraints that follow from differentiation with respect to time or from 
imposing gauge fixing conditions, but rather from adjoining the adjoint set of constraints. This is 
a feature of the Rarita-Schwinger constrained fermion system that has no analog in the familiar 
constrained boson systems such as gauge fields. 

Evaluating the brackets shows that M has the general form 

/0 -1 0 0 ^ 

1 U S T 

M = 

0 V A B 
\ 0 W C V j 

where in the SU(n) gauge field case, each entry in M is a 2 n x 2 n matrix (corresponding to the 
fact that 4> a is implicitly a 2 n component column vector, and Xb is implicitly a 2 n component row 
vector). Evaluating det M by a cofactor expansion with respect to the elements of the two unit 
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matrices ±1, we see that the submatrices U, S , T, V, W do not contribute, and we have 


detM =detJV 



(60) 


So we need to only evaluate the brackets M 33 = A, M 34 = B, M 43 = C, M 44 = V , giving 

A = — 2iga ■ B(x)5 3 (x — y) , 

B = -2Dg- Lg6 3 (x-y) , 

C =2Lg ■ D s 5 3 (x - y) , 

V =(i(L ?) 2 + B ■ (L s x L s ))5 z (x-y) . 

(61) 


When L = D, these become 


A = — 2 iga ■ B(x)5 3 (x — y) , 

B = -2{D s f5\x-y) , 

C=2(Dz) 2 5\x-y) , 

V =i((D ?) 2 - ga ■ B(x))fr i (x - y) . 

(62) 


Reflecting the fact that the and Xa ar e adjoints of one another, together with the fact that 
the matrix relating \Ed to P is anti-self-adjoint (see Eq. (1541) ) . these matrix elements obey the 
adjoint relations 

M ab (x,y) ] = ~M ba {y,x) . (63) 


Applications of these bracket and determinant calculations will be made in the subsequent paper, 
where we discuss quantization by both the Dirac bracket formalism and by the Feynman path 
integral. 

To conclude this section, we note that the constraints x, x' ■, TvT' 0 > -F^t play the role of gauge 
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transformation generators. For example, we have (with common time argument t suppressed) 


J d 3 x-x\x)e(x),$(y) 


- / d 3 xPq, o (x)D 0 e(x)^ 0 (y) 


= D -ije(y) , 


C 


C 


=Doyt(y) ■ 


(64) 


So the fermionic gauge transformation is a canonical transformation. This is also evident from the 
fact that since Eq. (13611 is just a shift in the fermionic variables T and To by the quantities De and 
D 0 e, which have no dependence on the fermionic variables, this shift leaves the canonical brackets 
[^i,Pj] c , [T 0 , P<£ 0 ] c etc. unchanged. 


V. GENERALIZED GAUGE INVARIANCE OF THE RARITA-SCHWINGER ACTION 

We turn now to a justification of our claim that the fermionic gauge transformation introduced 
in Eqs. ([9]) and (1361) is a generalized form of gauge invariance, which corresponds to redundant 
degrees of freedom, and which leaves essential attributes of the physics of gauged Rarita-Schwinger 
fields invariant. In the most familiar gauge invariant theories, such as Abelian or non-Abelian gauge 
fields, the Lagrangian density is invariant under a gauge transformation on the fields. These theories 
exhibit what one could term “strong” gauge invariance. In a weaker form of gauge invariance, which 
occurs for the free Rarita-Schwinger equation, the Lagrangian density changes by a total derivative 
under a gauge transformation of the fields, and so only the action is gauge invariant. Characteristic 
features of this case have been studied by Das [l2]. We argue in this section that there is a still 
weaker form of gauge invariance, obeyed by the massless Rarita-Schwinger equation with Abelian 
or non-Abelian gauging, in which under a gauge transformation the Lagrangian changes by a total 
derivative plus terms which vanish when initial value constraints are obeyed. 

We divide our argument that the transformation of Eqs. m and (1361) is a generalized form a 
gauge invariance into two parts, first considering infinitesimal transformations, and then considering 
general finite transformations. 


A. Infinitesimal gauge transformations 


In his seminal analysis of constrained systems, Dirac [141 ] classifies as “first class” constraints 
the maximal set of constraints that have vanishing mutual Poisson brackets, and notes that “Each 
of them thus leads to an arbitrary function of the time in the general solution of the equations of 
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motion with given initial conditions”. Elaborating on this, he notes that “Different solutions of the 
equations of motion, obtained by different choices of the arbitrary functions of the time with given 
initial conditions, should be looked upon as all corresponding to the same physical state of motion, 
described in various way (sic) by different choices of some mathematical variables that are not of 
physical significance (e.g. by different choices of the gauge in electrodynamics or of the co-ordinate 
system in a relativistic theory.)” 

These remarks suggest that gauge invariance, in its most general form, corresponds to an ar¬ 
bitrariness in the time evolution of a system, in the sense that the future evolution of the system 
is not uniquely determined by the initial conditions and the Euler-Lagrange equations following 
from the action principle. Under this generalized definition, the Rarita-Schwinger equation with 
coupling to gauge fields has a fermionic gauge invariance. To see this, we note the Euler-Lagrange 
equations yield equations of two types. The first are the time evolution equations contained in Eq. 
([3]). that determine the field variables at a later time t + At from those initially given at time t. The 
second are the primary and secondary constraints of Eqs. Q and Q, which constrain the initial 
field values at time t. If we make the gauge transformation of Eq. ([9|) at time t, with infinitesimal 
gauge parameter e (with its adjoint), we have seen that the action at time t changes, to first order 
in e, according to Eq. m- So assuming that the initial data at time t obeys both the primary and 
secondary constraints, then when the constraints at time t are applied the change in the action is 
0((e) 2 ). After this gauge transformation, we have seen in Eq. (1391) that the Euler-Lagrange equa¬ 
tions V, the primary constraint x an d the secondary constraint uj are all changed at order e, but 
because the gauge transformation is a canonical transformation that preserves inner properties, we 
have also seen that the altered secondary constraint is the one implied by the altered V and X- with 
an error of at most (e) 2 . Hence after the gauge transformation, we still have consistent equations 
of motion and initial conditions, which can serve as a starting point for time evolution. However, 
by making the gauge infinitesimal gauge transformation, we have introduced an arbitrariness into 
the evolved solution. In order to get a unique time evolution path from the initial data at time t 
using the action principle, one must impose a gauge fixing condition, that selects one member out 
of the equivalence class of equal action field configurations. 

In the gauged Rarita-Schwinger theory only the constrained action and constrained fermion 
number, in both flat and curved spacetimes, are invariant to first order under infinitesimal fermionic 
gauge transformations. This has an important physical significance. Consider a set of Rarita- 
Schwinger fields that, as envisaged in the model of Q], are permanently bound into condensates. 
The only way to see that these fields are present is through their gravitational fields, through their 
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gauge field polarizabilities, and possibly also through their influence on overall fermion number 
counting. The constrained action is the functional of the metric and the gauge fields that determines 
the influence of the Rarita-Schwinger fields on the metric and the gauge fields respectively, so the 
fact that the constrained action is invariant under infinitesimal fermionic gauge transformations 
means that the physical effects induced by confined Rarita-Schwinger fields are similarly invariant. 
(This statement is not contradicted by the fermionic gauge non-invariance of the energy integral 
and the gauge field source currents, since these are calculated by varying the unconstrained action, 
and do not take into account the fact that the constraints that enter into the constrained action 
are themselves non-trivial functions of the spacetime metric and the gauge fields.) 

The fermionic gauge invariance of the constrained action functional of the metric and the gauge 
fields then allows us to impose a gauge fixing constraint, making the time evolution determined by 
the action principle unique. Gauge fixing eliminates the redundancy of gauge degrees of freedom, 
and so is a convenience in checking the correct helicity counting for the Rarita-Schwinger fields, 
but is not needed for this purpose. In the following paper, where we turn to quantization, gauge 
fixing is needed to get an invertible constraint matrix in the weak field limit, and when covariant 
radiation gauge fixing is used one finds manifestly positive semi-definite anticommutation relations 
for the quantized Rarita-Schwinger fields. 

B. Finite gauge transformations: auxiliary fields and the extended action 

Since the transformations of Eqs. Q and (1361) are linear, and since the Euler-Lagrange equations 
and primary and secondary constraints are linear in the Rarita-Schwinger field, the relations of 
Eq. (|39D give the most general form of the transformed equations of motion and constraints. 
Thus, letting A denote a finite fermionic gauge transformation, the general form of the equations 
of motion and constraints are 

0 = 1?(A) =(? x o + D x T — a x -DoVk — ig{B + a x E) A , 

0 = x(A) =(? • D x T — iga ■ BA , 

0 = w(A) =<r ■ R(T 0 + DqA) - (B + a x E) ■ (T + DA) . 

( 65 ) 
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Under the gauge shifts of Eq. (1361) . A is augmented to A + e, or equivalently, under the extended 
gauge transformation that includes a shift of A, 

d'o ->■ d'o + D 0 e , T —» T + De , A A - e , ( 66 ) 

the formulas of Eq. (|65D are left invariant. By using Eq. (14311 . one can verify that 

D • V(A) - D oX (A) = igu{A) . (67) 

From Eqs. (1651) one deduces alternative forms of the \E r equation of motion, subject to the constraint 

X(A) = 0, 

£> 0 T =15^0 +iD x T + g(B - iE)A , 

0 = 6 ( A) = a ■ D x hp — Dqo • 'h — iga ■ EA . ( 68 ) 

From the first of these one finds 

D 0 D- V = {Df^o + g{B + iE) ■ $ + gD ■ ((B-iE)A) , (69) 

which gives a condition on the gauge shift A for the covariant radiation gauge condition D ■ 'h = 0 
to be maintained in time. 

We can now write down an action corresponding to the generalized equations of motion and 
constraints. It is 

S(A) J d 4 x[—'hjcr • D X ® + & ■ (a x D^ 0 + D x # - a x D 0 $) 

—igffi ■ (B 4 - <7 x E) ■ A + igA^(B + a x E) ■ 'h 

+ig x hgU • BA — igA^a ■ B'&o 

+igA^ (B + a x E) ■ DA — igA^a ■ BDqA] 

(70) 

One can check that the final line of this action is self-adjoint, by using Eq. (1431) . and one can 
also verify that this action is exactly invariant under the transformation of Eq. (1661) . including 
quadratic terms in e, without using the constraints following from the equations of motion. The 
extended action of Eq. (1701) gives the most general form of the gauged Rarita-Schwinger action, in 
which A plays the role of an auxiliary field that restores exact fermionic gauge invariance. 

Varying this action with respect to gives the generalized equation of motion U(A) = 0, while 
varying it with respect to Tq gives the generalized primary constraint x(A) = 0. Since these hold 
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for all times, Eq. (1671) then shows that they imply the generalized secondary constraint <u(A) = 0. 
Varying this action with respect to At gives just the secondary constraint oj( A) = 0 as the equation 
of motion for A. This shows that A is not an independent dynamical variable but rather is a 
Lagrange multiplier for the secondary constraint, which plays the role of a generalized auxiliary 
field. This further supports our argument that the gauge transformation of Eq. (1361) corresponds 
to a generalized gauge invariance, and that the gauge degrees of freedom are redundant degrees of 
freedom. 

Making the shift e = —A reduces A to zero, so that action of Eq. dZQD reduces to its first line, 
which is the original action of Eq. (1231) . Conversely, this shows that Eq. (j70l) is just Eq. (]23l) with 
the substitutions To — >• To + DqA and T — » T + DA, that is 

5(A) J d 4 x[-(Tj, + At ■ D x (T + DA) 

+(Tt + A*Id) • (a x D(To + D 0 A) + D x (T + DA) -ax D 0 (T + DA))] , 

(71) 


which makes manifest the invariance of S( A) under the shift transformation of Eq. (1661) . The 
simplicity of this way of constructing the extended action is a reflection of the fact that the 
fermionic gauge group is simply an Abelian group under addition of gauge functions. If we now 
define Tq = To + Do A and T 7 = T + DA, and fix the choice of A by imposing a gauge fixing 
condition, such as gauge covariant radiation gauge, then we see that as function of the primed, 
gauge-fixed variables the generalized action 5(A) takes the same form as the original action of Eq. 
(1231) took as a function of the original variables. 

The above analysis in terms of two-component, left chiral spinors can also be carried out in 
the original four-component formalism. Making the substitution t ^ + D /X A in Eq. Jl]) gives 
after some algebra using Eq. fl7]) the four-component form of the extended action functional of the 
Rarita-Schwinger field i/> p and the auxiliary field A, 


S{A) = l - J ^ 1 t]Dv^ p 

A —( — A'y^YiF^'tpp + 'i/’ /i 757r?D 1 /pA — A'y^'y r jF U pD p A)] 


(72) 


which is self adjoint by virtue of the Bianchi identity 


e^P[D p ,F vp \ = 0 


( 73 ) 
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Varying Eq. (1721) with respect to gives the generalized Euler-Lagrange equations (which include 
the generalized primary constraint) 

e^{D^ p + 9 -F up E) = 0 , (74) 

while applying g~ l D p to this and using Eq. d73l) gives the generalized secondary constraint 

e^ up F p M P + D p A) = 0 . (75) 

Varying Eq. (1721) with respect to A gives just the generalized secondary constraint of Eq. ()75l) , again 
showing that A is a Lagrange multiplier for the secondary constraint which acts as an auxiliary 
field, and thus corresponds to a redundant degree of freedom, not a physical degree of freedom. 


VI. PROPAGATION OF A RARITA-SCHWINGER FIELD IN AN EXTERNAL 
ABELIAN GAUGE FIELD: ABSENCE OF SUPERLUMINAL PROPAGATION 


We specialize now to the case of a Rarita-Schwinger spinor propagating in an external Abelian 


gauge field, as studied by Velo and Zwanziger [4], For an Abelian gauge field, 

1 <7 ■ B 


(76) 


a-B (B) 2 

and so a ■ B is invertible as long as ( B ) 2 > 0, which we assume. Provided the Lorentz invariant 
expression ( B ) 2 — ( E ) 2 is positive, ( B ) 2 will be positive in any Lorentz frame. In discussing 
undamped wave propagation we will not use the inequality ( B ) 2 — ( E) > 0 , but in treating damped 
longitudinal mode propagation in Appendix B, we will assume that (E ) 2 /(B ) 2 is small, as motivated 
by the fact that when ( E ) 2 is of order (B ) 2 the vacuum is highly unstable against pair creation. 
(Strictly speaking, the vacuum is stable against pair production only when E ■ B = 0 and ( B) 2 — 
(E ) 2 > 0, that is, when there is a Lorentz frame in which the Abelian field has vanishing E 115].) 
Given that (B ) 2 > 0, we can solve the constraint u = 0 of Eq. (1251) for 4'o, giving 

Q-$ 


4 / n = 


(BY 


(77) 


where we have defined 


Q = a ■ B(B J t~axE) = BxE + Be j • (B + iE ) — iB ■ Ea . (78) 

Substituting the solution for 'Lq into Eq. (13TT) . we get an equation of motion for 'L by itself, 

D 0 ^ = D-^— + iDx 'L . 

(B ) 2 


( 79 ) 
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To determine the wave propagation velocity in the neighborhood of a spacetime point x* = 
(i*,x*), we need to calculate the equation for the wavefronts, or characteristics, at that point. 
Writing the first order Eq. (179j) in the form 

<9 0 T = v^: 

( 5 * 

with B * and Q* the values of the respective quantities at x*, we see that A[\I/,x*,x] involves no 
first derivatives of T at x*, and so is not needed [1(|, [17] for determining the wavefronts of Eq. 
(1311) . The reason is that when taking an infinitesimal line integral of Eq. (1801) . according to 

lim f d£[d oT = ...] , (81) 

J-6 

discontinuities across wavefronts contribute through the first derivative terms, but when the exter¬ 
nal fields are smooth the term A[T,x*,x] makes a vanishing contribution as 5 —>■ 0. Dropping A, 
and multiplying through by (-B*) 2 , we get the equation determining the wavefronts in the form 

(B*) 2 d 0 $ = VQ* • T + f(i?*) 2 V x T . (82) 

By similar reasoning, the constraint x can be simplified, for purposes of determining the wavefronts, 
by replacing T by V, giving 

0 = S ■ V x T . (83) 

Since these are now linear equations with constant coefficients, the solutions are plane waves, 
and without loss of generality we can take the negative z = x 3 axis as the direction of wave 
propagation. So making the Ansatz 

T = (7exp(fDf + i Kz) , (84) 

Eq. (1821) for the wavefronts or characteristics takes the form 

0 = F= {B*) 2 nC - KzQ * • C - i(B*) 2 Kz x C , (85) 

with z a unit vector along the z axis, and the constraint Eq. (1831) becomes an admissability 
condition on C, 




+ iV x f + A[\k, x*, x] 


(80) 


0 = a ■ z x C 


(86) 
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Writing F m as a matrix times C n (and dropping the subscripts *, which are implicit from here 
on) we have 


Fm, — -\ r rn n Cn 


Nmn — (B) C^inn F5 m 3Qn 


i(B') Ke m 3 n 


(87) 


The equation for the characteristics is now 


det(iV) = 0 , ( 88 ) 

since this is the condition for Eq. f(85f) to have a solution with nonzero C. However, since evalua¬ 
tion of the determinant shows that it factorizes into blocks that determine C\ 2 and a block that 
determines C 3 , a simpler way to proceed is to work directly from the equations F m = 0 , which 
decouple in a corresponding way. Calculating from Eq. 11851) . we find 

0 =Fl’ 1 = (B) 2 (ftC'J ' 4 + iKC\ A ) , 

0 =F^ 1 = (B) 2 (yLCl A - iKC\ A ) , 

0 =f } 1 = (B) 2 nc} 1 -K(Q-C)^ 1 , 

(89) 

where "f, -I indicate the up and down spinor components, labeled in Eq. (I22[) by a = 1, 2. Similarly, 
the constraint Eq. () 86 l) becomes 0 = — CT 1 C 2 + 02 C 1 , that is 

C$=iCl , 

ni — __ Ani 
L / 2 — , 

(90) 

with no corresponding condition on C([ : ^. The first two lines of Eq. (1891) together with Eq. (1901) 
have the solution 

cl =C , cl = iC , Q = K , 
cl =C , C'l = -iC , n = -K , 

(91) 


with C arbitrary, corresponding to waves with velocity of magnitude |fl/A'| = 1. Thus the modes 
with Ci , 2 7 ^ 0 are exactly luminal. Because general background gauge fields are a non-isotropic 
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medium, these modes have nonzero longitudinal components given by solving the third line of Eq. 

(USD, 

C 3 = K((B) 2 n-KQ 3 )-\Q 1 C 1 + Q 2 C 2 ) . ( 92 ) 

The effect on the characteristics of a gauge change T — > ^ + De , e = E exp(ifR + iKz)f(t, z ), 
where / has a unit slope discontinuity along the z axis at x*, is to shift —> C'J’'*' + and 

thus are gauge degrees of freedom. In Appendix B, we continue this discussion and show that 
the longitudinal gauge mode with C\ = C 2 = 0, C 3 ^ 0 also does not propagate superluminally, 
although in general it is subluminal. 


VII. FAILURE OF ADIABATIC DECOUPLING AND INAPPLICABILITY OF THE 

5-MATRIX “NO-GO” THEOREMS 


We show in this section that various “no-go” theorems that claim to rule out gauging of higher 
spin theories do not apply to the gauged Rarita-Schwinger field. The reason is that there is a failure 
of adiabatic decoupling, arising from the fact that the uj secondary constraint is homogeneous in 
the gauge fields. For a recent paper on “no-go” theorems see [lO], which has extensive references 
to the earlier literature. In our analysis here we shall refer specifically to the paper of Porrati jj|, 
which uses so called “on-shell” methods to give limits on massless high-spin particles. 

The analysis of Porrati assumes that “the general helicity-conserving matrix element of a 17(1) 
current between on-shell spin s states is {v,p + q\J fl \u,p )... v , where u and v are free-space spinors 
that obey the massless Dirac equation. Porrati assumes that the matrix element is bilinear in u 
and v , and “otherwise depends only on the momenta”. We shall see in the following subsections 
that this assumed form is not realized in the gauged Rarita-Schwinger theory, where because of the 
failure of adiabatic decoupling the matrix element in question also depends on the 17(1) gauge field 
polarization through the dual field-strength F /w = ^e^x^F^. In fact, the initial and final Rarita- 
Schwinger spinors both must have a 7) x „ dependence in order to obey the secondary constraint of 
Eq. Q, and so the matrix element has the more complicated form (v,p + q , F^l J^\u,p, F^ v ). 

We show in Sec. 7A that the initial and final Rarita-Schwinger spinors in the limit of zero gauge 
field amplitude are equal to free-space spinors u, v of the form assumed by Porrati, plus a fermionic 
gauge transformation that depends explicitly on the photon field strength F^ v . This structure 
arises from the homogeneous form of the secondary constraint, and corresponds to an intrinsically 
non-perturbative aspect of the gauged Rarita-Schwinger equation. As another reflection of this, 
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we show in Sec. 7B that one cannot set up a covariant Lippmann-Schwinger equation [l8] for the 
Rarita-Scwhinger wave function, and so the matrix element that enters into the “no-go” theorems 
does not admit a Born approximation. In Sec. 7C, we show that a matrix element that has all the 
required invariances can be formulated using an analog of the distorted wave Born approximation, 
in which the initial and final Rarita-Schwinger states have an explicit dependence on the photon 
polarizations. 

A. The zero amplitude limit of the \F equation: retained memory of the gauge field 

As in Sec. 6, let us consider a Rarita-Schwinger field propagating in an external Abelian gauge 
field. For convenience, we assume that the ratio \E{x)\/\B{x)\ = r{x) is bounded from above. In 
the limit as the vector potential amplitude A is scaled to zero, Eqs. (1771) and (178]) become 

'F 0 (£) =R{x) ■ $(x) , 

R(x) =a ■ B(x)(B(x) + r(x)a x E(x)) , 

(93) 

with B = B/\B\ and E = E/\E\ unit vectors along the E and B fields. When the external field 
is a propagating plane wave with wave vector direction q , the unit vectors q , B and E form an 
orthonormal set of constant unit vectors, and |r(x)| = 1. We see that because the secondary 
constraint of Eq. ([6]) is homogeneous in the field strengths, the relation between 'Fg and \F retains 
a memory of the gauge field orientations, and thus of the photon polarization, even in the limit as 
the field amplitude approaches zero. 

In the zero amplitude limit, Dq = do and D = V, so substituting Eq. (1931) into Eq. (1791) . the 
zero amplitude limit for the equation of motion for \F becomes 

do'F = VR -f-MVxf . (94) 

with the primary constraint now a • V x *F = 0. Hence through R the ^ equation of motion retains 
a memory of the external fields in the limit of zero amplitude, that is, adiabatic decoupling has 
failed. Let us now consider the situation in which the Rarita-Schwinger field and the external 
gauge fields are plane waves, so that R is a constant and \F has the form 


^ _ (j e i(nt+k-x) 


( 95 ) 


28 


Making the fermionic gauge transformation 

T T' = T + Ve , 

e = E e i(Clt+k-x) 

( 96 ) 

IL' still obeys the zero amplitude primary constraint since a ■ V x Ve = 0. Then the gauge choice 

E = (97) 

R ■ k 

reduces Eq. (IM1) to the free-space form 

<9 0 T' = iVx$' . (98) 

Thus a Rarita-Schwinger plane wave in a zero amplitude gauge field plane wave background is 
equal to a free-space solution plus a gauge term that has a memory of the photon polarizations. 

B. Breakdown of the Lippmann-Schwinger equation: no Born approximation to scattering 

Let us now examine what happens if one tries to set up a covariant Lippmann-Schwinger 
equation, so as to generate a Born perturbation series for the Rarita-Schwinger wave function 
in an external gauge field. Let us start from the Rarita-Schwinger equation in the form (see Eq. 

(USD) 

= 0 . (99) 

Splitting D v into d v and gA u , this equation takes the form 

i vvp d^ p = - 7 vup gA^ p . ( 100 ) 

Let us now try to solve this equation as a perturbation series around a free-space solution by 
writing 

^p(®) = ^p ee ( x ) + j d t yS pa (x-y)'f pK gAp(y)il; K (y) , (101) 

where ^ ree obeys the free-space Rarita-Schwinger equation 

y’^Vf 66 = 0. (102) 

If the free-space Green’s Rarita-Schwinger Green’s function S pa (x — y) obeyed 


'y vup d x „S pa (x - y) = -5 p 5 4 (x - y) , 


( 103 ) 
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then Eq. (11011) would reproduce Eq. (HOOD . But in fact the free-space Green’s function cannot 
obey Eq. (11031) . because d xri r f lvp d xv S p(x [x — y) = 0; instead it obeys [221 ] 


'f' vp d xv S pa (x - y) = -5%5 a (x - y) + d ya Q v (x - y) 


(104) 


with O necessarily nonvanishing. Integrating d ya by parts onto the factor r y apK gAs{y)ip K (y), 
gets 


one 


7 al3K gF a g{y)'ip K (y) + ^ K gAp(y)d ya ^ n {y) 


(105) 


The first term of this expression vanishes by virtue of the secondary constraint, but the second 
term is non-vanishing because the Rarita-Schwinger equation for the exact wave function ip K {y) is 

^ K D m My) = o , (foe) 

that is, it requires the full covariant derivative D ya in place of its free-space restriction d ya . The 
conclusion from this analysis is that one cannot set up a covariant Lippmann-Schwinger equation 
for the gauged Rarita-Schwinger wave function, and thus one cannot develop this wave function 
into a Born approximation series expansion in powers of the coupling g to the external gauge field. 


C. Lorentz covariance and mode counting in on-shell Rarita-Schwinger field-photon 
scattering: a distorted wave Born approximation analog 

We address finally the question jll]] of whether one can write down an amplitude for leading 
order on-shell scattering of Rarita-Schwinger fields from an external electromagnetic field, which 
has the requisite relativistic covariance while preserving the correct counting of massless spin | 
propagation modes. Looking ahead to quantization, an operator effective action for this scattering 
process can be inferred from the interaction term in Eq. ©, 

Ses{^Pfii j 4j/) — J d x(L e fi ('i/’/x, A v ~) , 

AffW’/oA') =^g^ p {x)ie mup ^ v A u (x)iip(x) , 

(107) 

where we have suppressed spinor indices as in the text from Eq. (]3]) onwards. For Abelian external 
fields A u , the covariant derivatives in the equations of motion and constraints are given by 

D u = d u + gA u , ^5 U = < d u - g A u 


(108) 
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At the outset we shall assume that A u (x) is of short range, and vanishes for \x\ > R for some 
radius R. This effective action, the equations of motion of Eqs. Q and (|4|), and the primary 
and secondary constraints following from them, given in Eqs. (J5]) and (|6]), are all relativistically 
covariant, and so provide a starting point for calculating a covariant scattering amplitude. Taking 
the matrix element of Eq. (110711 between an incoming Rarita-Schwinger state of four-momentum p, 
and an outgoing Rarita-Schwinger state of four momentum p' , we get the corresponding scattering 
amplitude 

• As = \ ia j ’ x ) e ^ P ^v A v(x)i> p {Pix) , (109) 

where i^p and are now wave functions, rather than operators, that obey the Rarita-Schwinger 
equations of motion in the presence of the external held A u . 

We now introduce source currents for the gauge potential A u and the Rarita-Schwinger wave 
functions ip p and ip p , and study their conservation properties. The source current to which the 
gauge potential A u couples is defined by writing the scattering amplitude as 

As =\ig J d 4 xA u (x)J u (x) , 

J v (x) ='ip ll (p',x)e firiup 'y 5 'y ri 4’ p (p,x) . 

( 110 ) 


The source current for the Rarita-Schwinger held ipfi(p'i x ) is dehned by writing the scattering 
amplitude as 


As=-ig J d i xtp p {p\x)J p {p,x) 
J P {p,x) =e /W 757 r,A u (x)4’p(p,x) . 


(Ill) 


Finally, the source current for the Rarita-Schwinger held V’p(P) x) is dehned by writing the scattering 
amplitude as 




^ig J d 4 xJ P (p, x)il) p {p, x) 


J P {p',x) =i4 p {p' ,x)e pr]Vp ^ r] A v [x) 


( 112 ) 


We now show that the three currents that we have just dehned are conserved. For the source 
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current J v for the gauge potential, we have 

duJ u =ipp{p',x)5 v d ir)vp ^ ri ip p {p,x) 

+^p{p, x)e wl/p 757 r,D v ip p {p, x) 

=0 , 

(113) 

where the first and second terms on the right vanish by the Rarita-Schwinger equations for ip „pp', x) 
and ip p (p,x) respectively. For the source current J p {p,x) for the spinor ip (p',x) , we have 

DpJ p {p, x) =e WI/p 757 r, (dpA u (x))ipp(p, x) 

+e w/p 7 5 7 r} A v {x)D ll ip p {p, x) 

=0 . 

(114) 

The second term on the right vanishes by the Rarita-Schwinger equation for ip p (p,x), while the 
first term on the right can be rewritten as 

^e wl/p 757 v Fp U (x)ipp(p,x) (115) 

and vanishes by the secondary constraint of Eq. ([6]). Finally, for the source current J P (p', x) for 
the spinor ip p (p,x), we have 

J p (p', x)J5p =ipp(p', x)e^ p 7 5 7 v (d p A v (x)) 

+ipp,{p', x)l5 P e pr]Vp ^r ] A v [x) 

=0 . 

(116) 

Again, the second term on the right vanishes by the Rarita-Schwinger equation, while the first 
term on the right vanishes by the secondary constraint of Eq. d6j) . 

Consider now the following three gauge transformations, 

A v {x) — >A v {x) + d u A , 
fp P (p,x) -+ip p {p,x) + D p a , 
ipn(p', x) -^ipp(p, x) + pbp, , 


(117) 
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with a and f3 independent spinorial gauge parameters. From Eqs. (IIlOD - dl 121) . together with Eqs. 
(I113I) - (I116I) . we find that these transformations each leave the amplitude As invariant, 

=\ig j d 4 x[d u A)j u (x) = — ^9 j d 4 xAd v J u (x) = 0 , 

d a As =\ig J d 4 xJ p (p\x)D p a = ~ig J d 4 xJ P (p', x)*D p a = 0 , 

bpA s =\ig J d 4 xfj i D ll J ll (p, x) = ig J d 4 x(3D p J p (p,x) = 0 . 

(118) 

This, together with the primary and secondary constraints, implies the correct mode-counting 
for the Rarita-Schwinger wave functions, since the gauge degrees of freedom do not change the 
amplitude and so are redundant. 

We next must specify more precisely the structure of the spinor wave functions entering the 
formula for As- Since the gauge field A v is assumed to vanish in the external region |x| > R, the 
Rarita-Schwinger wave functions obey free field equations in this region. So for \x\ » R they can 
be taken asymptotically as plane waves at t —>• ±oo, 

i> l i(p',x) ~u fl (p')e ip '' x , t } Too , 
ip P (p,x) ~ Up (p)e ip ' x , i->-oo . 

(119) 

With these boundary conditions, the formula for the amplitude takes the final form 

As = \ ig J d 4 ^ \p'i x )t MUp lblr 1 A v{ x )^\p, x ) ■ (120) 

The out state (-) and in state (+) boundary conditions used here are analogs of the boundary 
conditions used in the distorted wave Born approximation [20], which the construction of Eq. (11201) 
resembles. Equation (11201) then gives an approximation to the matrix element for Rarita-Schwinger 
scattering by the gauge potential. 

Rather than invoking the presence of redundant degrees of freedom to count physical Rarita- 
Schwinger states, we can follow the usual procedure of imposing a gauge-fixing constraint. To 
preserve relativistic and gauge covariance, this can be taken as the gauge covariant Lorentz gauge 
condition 

®)^ M = D p ^ p (p, x) = 0 . (121) 

which is attainable from a generic gauge by the gauge transformation of Eq. m, provided that 
D p Dp is invertible. In the external region where the gauge field vanishes, one can instead use 













33 


the condition 7 p i\) p = 0 in place of the secondary constraint together with the gauge condition 
d p i\) p = 0 , giving the usual covariant degree of freedom counting for the incoming and outgoing 


Rarita-Schwinger wave functions 


21] - Alternatively, if we are not concerned to maintain manifest 


Lorentz covariance, we can make a gauge transformation in the external region to the gauge 
V’o = V • = 0 used in 


mi, H inger degrees of freedom. When a 

non-Lorentz covariant radiation gauge condition is used, scattering matrix elements depend on a 
unit timelike vector in addition to the particle momenta, and so the conditions assumed in {{J are 
not obeyed. 

Note that if one were to attempt to construct a Born approximation amplitude, in which the 
Rarita-Schwinger wave functions in the presence of the gauge field are replaced by plane waves in 
the interior region where the potential is nonzero, the arguments given above for compatibility of 
Lorentz covariance with degree of freedom counting would fail. The reason for this is that the spinor 
source currents would then no longer be conserved, even to zeroth order in the gauge coupling g, 
because the free particle plane wave solutions do not obey the secondary constraint of Eq. Q. The 
non-existence of a satisfactory Born approximation for Rarita-Schwinger photon scattering agrees 
with the result obtained in Sec. 7B, that one cannot construct a Lippmann-Schwinger equation 
for this process. To establish compatibility, we have had to use an analog of the distorted wave 


Born approximation [201 ], in which the leading approximation to the amplitude is constructed using 
interacting rather than free fermion wave functions and does not have a perturbation expansion 
for small coupling, g . 

When the external Abelian potential is a plane wave field which extends to infinity, there is no 
large \x\ region where the Rarita-Schwinger solutions reduce to free-space ones. Rather, as shown 
in Sec. 7A, in the adiabatic decoupling limit of a zero amplitude gauge field, the Rarita-Schwinger 
solutions become free-space solutions plus gauge terms that remember the photon polarization, and 
which are necessary to enforce the secondary constraint. Thus one cannot attain the kinematic form 
assumed in the on-shell “no-go” theorems. But as shown here, using distorted Born approximation 
waves one can write down a consistent covariant scattering amplitude. 


VIII. SUMMARY AND REMARKS 

To conclude, we see that unlike the massive case, the massless gauged Rarita-Schwinger equation 
leads to a consistent classical theory. The theory has the correct counting of propagating non-gauge 
degrees of freedom with no superluminal wave propagation. The theory admits a generalized 
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fermionic gauge transformation, and infinitesimal gauge transformations are an invariance of the 
constrained flat and curved spacetime actions and of the fermion number. The gauged Rarita- 
Schwinger equation has a non-perturbative aspect when the secondary constraint u is eliminated, 
resulting in a breakdown of adiabatic decoupling, leading to the inapplicability of various S'-matrix 
“no-go”theorems that claim to forbid gauged massless Rarita-Schwinger fields. The extension of 
these results to the quantized Rarita-Schwinger theory is given in the following paper, where we 
show that a consistent quantization by the Dirac bracket and path integral methods is possible, 
with a manifestly positive semi-definite canonical anticommutator in covariant radiation gauge. 
Thus, in the massless case our analysis eliminates the various objections that have been raised 
to gauging Rarita-Schwinger fields, showing that non-Abelian gauging of Rarita-Schwinger fields 
can be contemplated as part of the anomaly cancelation mechanism in constructing grand unified 
models. 

We conclude with several remarks: 

1. We have introduced gauge fixing to make time evolution of the Rarita-Schwinger fields 
unique, but the analysis of this paper does not require gauge fixing. Specifically, if gauge 
fixing is not imposed, the correct helicity counting is still obtained because fermionic gauge 
degrees of freedom are redundant degrees of freedom, and are not physical. Gauge fixing 
makes this redundancy manifest by providing a condition that excludes the gauge degrees of 
freedom, but in analogy to the case of Maxwell electrodynamics, gauge fixing is not needed 
to get the correct physical state counting. On the other hand, in the following paper, where 
we turn to quantization, gauge fixing is needed. This can already be anticipated from the 
form of the constraint matrix N of Eq. ()60l) , which when gauge fixing is omitted reduces to 
the single element A = —2 iga ■ B(x)5 3 (x — y) which is not invertible in the small B limit. 
Inversion of the constraint matrix does not enter into the calculations of this paper, but is 
needed in the following paper both for Dirac bracket and path integral quantization. 

2. A possible exception to the non-perturbative behavior detailed in Sec. 7 is when the E and 
B gauge fields are random, since if Eq. (1771) is replaced by an average, denoted by AV, 

( * 0)AV ~ ( (J^)av ‘ ( ^ )AV ’ (122) 


it becomes 


(^o)av — ' Wav 


(123) 
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which is compatible with (TAav = • (T)av = 0, the customary free Rarita-Schwinger 

constraints employed in 11], [22]. This heuristic observation suggests that Rarita-Schwinger 
fields coupled to quantized gauge fields with zero background gauge held may have a per¬ 
turbative g —> 0 limit. 


3. In showing in the Abelian case that there is no superluminal propagation, the inversion of 
a ■ B to get To only required ( B ) 2 ^ 0. In the non-Abelian case, where B is itself a matrix, 
the conditions for invertibility are nontrivial and have yet to be analyzed. We will see in the 
following paper that this issue is side-stepped when the constraints are dealt with by the 
Dirac bracket or path integral procedures, since these do not require inversion of a ■ B when 
a gauge constraint is included. 
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Appendix A: Notational conventions and useful identities 


We follow in general the notational conventions of the book Supergravity by Freedman and 
Van Proeyen 22], The metric r]^ u is (—,+,+,+) and the Dirac gamma matrices 7 ^, 7 ^ obey the 


Clifford algebra 


(Al) 


Ifilu + 7i'7m — 
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They are given in terms of Pauli matrices 07 by 


7o = -7 = 


Ifj = l 3 = 


75 = *70717273 = 


0 -1 
1 0 


0 a 




°3 0 


1 0 
0 -1 


We also note that 


, _ ,0123 _ 1 

€0123 — — € — 1 , 


the left chiral projector Pi is given by 


-Pl =-(I + 75) , 


and the spinor ip is defined in terms of the adjoint spinor ip' by 

ip = ip'i'y 0 


(A2) 


(A3) 


(A4) 


(A5) 


As noted in [[22], the Rarita-Schwinger equation of motion can be written in a number of 
equivalent forms. When ordinary derivatives are replaced by gauge covariant derivatives, these are 
the vector-spinor equations 

e^ p lv D u iP p =0 , 
l wp D v 'P p =0 , 

7 p (D u ip p - D p ip v ) =0 , 

7 a D a {D a ip v - D u ip a ) =y([D p ,D rT \ip„ + [D u ,D p ]ip a + [D a , D u ]ip p ^j , 

(A 6 ) 

with only the fourth line, which is quadratic in the covariant derivative, involving more than just 
a substitution d v D u in the formulas of 22], Using r y v y up = 2y p , these also imply the spinor 


equation 7 up D u ip p = 0. These formulas play a role in verifying stress-energy tensor conservation, 


as does the identity 


23| 


0 — e' K<JIW {A T B\C rJ D p E y + A v B r C\D,jE p + A p B u C T D\E,j + A a B p C u D r E\ + AxB a C p D u E T 


(A7) 





37 


with A t , B\ , C a , DE u five arbitrary four vectors. This identity follows from 

0 = 8^e Xa ^ + 5 v T e aXatl + 5^e vaXa + 5"e^ aX + 5 x e a ^ a , (A8) 

which is easily verified by noting that A, <r, v must take distinct values from the set 0,1,2,3, 
and that r must be equal to one of these values. 

The fundamental identity for the Pauli matrices is 

^ah T itabc&c , (A9) 

with ei 23 = 1 and with the index c summed. We repeatedly use the following two identities that 
can be derived from Eq. (1A9I) . for a general three vector A that is proportional to a unit matrix 
in the spinor space and so commutes with 3, 

3 x (<r x A) = — 2 A + iff x A , 

(A x 3) x 3 = — 2 A + iA x 3 . 

(A10) 

Additional useful identities are 


3 x 3 =2i3 , 

33 - A =A — ia x A , 

3 ■ A 3 =A + ia x A , 

{3 x A) ■ 3 = — 2i3 ■ A , 

3 ■ (3 x A) =2i3 ■ A , 

&a&b = 2(<^afe 2 ®b&a) i 

B = iA — A x 3 hA = —(B x 3) 

(All) 

Gauge field covariant derivatives are 

Dp = d^, + gA^ , (A12) 

with the gauge potential A p = A^a and the gauge generators anti-self-adjoint, and with the 
components A') self-adjoint. The non-Abelian generators tA obey the compact Lie algebra 


[tA,tB\ = fABCtc ; 


(A13) 
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in the Abelian case we replace tA by — i- In writing field strengths E and B we pull out an 
additional factor of i to make them self-adjoint, so that we have the identities 


D x D = — igB , 

[D, Do] = - igE . 

(A14) 

We will also write a right-acting three-vector covariant derivative as = \^ + gA, and define a 
left-acting three-vector covariant derivative as ^ — gA, so that we have the integration by 

parts formulas 

J dPxAJ^^B = — J d 3 xAiDgB , 

DgS 3 (x - y) = - 5 3 {x - y)5$ . 

(A15) 


<— 


An analogous definition is used for the operators L and L which enter the gauge fixing condition. 

At the classical level, variables will be either Grassmann even or odd. Irrespective of the 
Grassmann parity of monomials A and B, the adjoint operation is defined by 22] 


(AB)t = 5+At 


For classical brackets, we follow the convention of Henneaux and Teitelboim 



/a FdG__dF^dG\ f d L Fd L G d L Fd L G \ 

\ dq l dpi dpi dq l J \ dO a di r Q 97r Q dQ a J ’ 


(A16) 


(A17) 


with €f the Grassmann parity of F, with d L a Grassmann derivative acting from the left, and 
with q i , pi ( 9 a , ir a ) canonical coordinates and momenta of even (odd) Grassmann parity. Using 
the classical bracket, the Dirac bracket is constructed from the constraints as discussed in Sec. 
2 of the following paper. To make the transition to quantum theory, the quantum commutator 
(anticommutator) is defined to be ih times the corresponding Dirac bracket (with h = 1 in our 
notation). Classical canonical brackets are always denoted, as above, by a subscript C, with a 
subscript D used for the corresponding Dirac bracket. We use the standard notations [A, B] = 
AB — BA for the commutator and {A, B} = AB + BA for the anticommutator. 
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To calculate the Dirac bracket, we use block inversion of a matrix. Let 

. Ai A 2 

M = 


M~ l = 



(A18) 


with Ai,...,A 4 themselves matrices. Then when A 4 is non-singular, the blocks Bi,...,B 4 of M 1 
are given by 


A 

=A\ 

— A 2 A 4 3 Az 

Bi 

=A‘ 

-1 

5 

b 2 

= - 

A 1 A 2 A 4 1 

Bs 

= - 

a^a 3 A- 1 

B 4 

=*k 

_1 + a 4 3 a 3 a 


(A19) 


Even though the blocks are noncommutative, Eqs. (1A18I) and (1A19D give an inverse that obeys 
M~ l M = MM' 1 = 1 . 

When the constraints (j) a and \a are combined into an 8 element set of constraints K a = 
0m Ka+4 = Xa, cl = 1,..., 4 then the bracket matrix S a b(x,y) = [K a (x), Hb(y)]c can be expressed in 
terms of the matrix M a h ( x , y ) of Eq. (l58l) as 


S(x,y) = 


0 M(x, y) 
M T (y,x) 0 


(A20) 


where Mj b (x,y) = Al^ a {x,y) is the matrix transpose. Defining the inverse M 1 (x,y) that obeys 
f d 3 zAf- 1 (x, z)M(z, y) = f d 3 zAI(x, z)M~ 3 (z, y) = 5 3 (x — y), it is easy to verify that 

0 M T ~ 1 (y,x) 

M -1 (x, y) 0 


S (x, y) = 


(A21) 


Appendix B: Analysis of the Rarita-Schwinger field in an external Abelian gauge field: 
propagation of the longitudinal gauge mode 


We continue here the analysis begun in Sec. 5 to study propagation of the longitudinal gauge 
mode. We must now solve for C'J’ ^ starting from Eq. (|89D with 2 = 0, so the third line of Eq. 






(1891) simplifies to 
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0 ={B) 2 ncl’ 1 - A"(g 3 C 3 ) t ’ 4 , 

Q 3 =B\E 2 — B 2 E\ + B 317 ■ (B + iE ) — iB ■ E 03 

(Bl) 


Writing this as 

M_/ u n E/12 \ l cA 
W Ve / 21 U22) \c 3 *J 

we find for the matrix elements 


(B2) 


U u =(B) 2 n - K[B X E 2 - B 2 E l - i^B^ + B 2 E 2 ) + B 2 ] , 

U 22 =(Bj 2 n - I<[B^ - B 2 E\ + i(B x E 1 + B 2 E 2 ) - B%\ , 

U 12 = — KBzlBi + iE\ — i(B 2 + iE 2 )\ , 

U 2 \ = — KB^lBi + iE\ + i(B 2 + iE 2 )\ 

(B3) 

The equation 0 = det(U) = U 11 U 22 — U 12 U 21 reduces, after dividing by an overall factor of ( B ) 2 , 
to 

0 = (B) 2 Q 2 - 2QK(B 1 E 2 - B 2 Ei) + K 2 (E\ + El-Bl) , (B4) 


with the solution 

n x± y 1 / 2 

K ~ (B )2 ’ 

X =B\E 2 — B 2 Ei , 

Y =(B!E 2 - B 2 E x f - ( Bf(E\ YE 2 - B I) . 

(B5) 

The analysis of the solutions of Eqs. (1B4D and (1B5[) divides into two cases, according to whether 
the roots of Eq. (lB5j) are both real, or both complex. The roots are both complex if 

{BiE 2 - B 2 E x ) 2 < ( B) 2 {E 2 + E 2 - B |) , (B6) 


which can be rearranged algebraically to the form 


[(B) 2 - (E 2 + E 2 )]B 2 < (B 2 + B 2 )(E 2 + E 2 ) cos 2 0 


(B7) 



where we have written 
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B x E 2 - B 2 E 1 ={B\ + Bl) l / 2 (El + E 2 ) 1 ' 2 sine/) , 

B\E\ + B 2 E 2 =(B 2 + Blf/ 2 {E\ + El) 1 ' 2 cos (j) . 

(B8) 

Since the right hand side of Eq. ()B7I) is non-negative, when the left hand side is negative the 
inequality is satisfied, and both roots are complex. Hence a necessary (but not sufficient) condition 
for both roots to be real is 

(B) 2 - (El + El)> 0 . (B9) 

1. The hyperbolic case: both roots real 

When both roots are real, Eq. (1B1I) describes the hyperbolic case of propagating waves. Intro¬ 
ducing the velocity V = £l/K, Eq. (1B4I) can be written as 

0 = (B) 2 V 2 -2V(B 1 E 2 -B 2 E 1 ) + El + El-Bl , (BIO) 

which can be rearranged algebraically to the form 

[(Bl + Blf /2 -(El + Elf /2 } 2 + (B) 2 (V 2 - 1) = 2(Bf + Blf /2 (Ef + Elf /2 (Vsin . (Bll) 

Let us now assume that V 2 > 1, and show that this leads to a contradiction. When V 2 > 1, 
the left hand side of Eq. (IB111) is nonnegative, which implies that V sin 0 on the right must be 
nonnegative, and so can be replaced by its absolute value. Hence the right hand side of Eq. (IB111) 
obeys the inequality 

2(Bl+Blf/ 2 (El+Elf /2 (Vsincl)-l) = 2(H 1 2 +H|) 1/2 (Ef+E|) 1/2 (|Hsin</.|-l) < 2(H) 2 (|H|-1) , 

(B12) 

where we have used Eq. (IB9D . But the left hand side of Eq. (1B11D obeys the inequality 

[(Bl + Blf/ 2 - (El + Elf/ 2 } 2 + (Bf(V 2 - 1) > (H) 2 (|H| + 1)(|H|-1) > 2(H) 2 (|H|-1) , (B13) 

which is a contradiction, since a real number cannot be strictly less than itself. Hence we must 
have V 2 < 1, and there is no superluminal propagation. 
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2. The elliptic case: both roots complex 


When both roots are complex, Eq. m describes the elliptic case in which there are no 
propagating waves; when a propagating wave enters an elliptic region from a hyperbolic one it will 
be damped to zero amplitude. However, in the case of weak damping one can still define a wave 
velocity and ask what its magnitude is. When both roots are imaginary, Eq. (1B5I) takes the form 

n _x±i(-Y) 1 / 2 3 4 5 6 7 8 

k~ > 

X =B\E2 — B 2 E 1 , 

—Y = - (BtE 2 - B 2 E 1 ) 2 + (B) 2 (E 2 + El - B 2 ) . 

(B14) 


Regarding H as real and the wave number K as complex, the effective propagation velocity has 
the magnitude 




n 

~x~r 


X 2 - Y _ E\ + E%-Bl 
(B) 2 \X\ ~ \B 1 E 2 -B 2 E 1 \ 


The condition for weak damping is — Y « X 2 , which can be rewritten as 


(B15) 


(B) 2 (E% + E 2 - B 2 ) « 2{B 1 E 2 - B 2 E 1 ) 2 , (B16) 


and implies 


|Reff| « 


2\B\E2 — B 2 E 1 

w 


2| E 
\B\ 


(B17) 


Hence as long as 2\E\ is not much larger than \B\, which is required by the vacuum stability 
condition \E\ < \B\, the damped wave propagation velocity is subluminal. 
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